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We analyze the spectral energy distributions (SEDs) of Lyman break galaxies 
(LBGs) at z ~ 1-3 selected using the Hubble Space Telescope ( HST ) Wide Field 
Camera 3 (WFC3) UVIS channel filters. These HST /WFC3 observations cover 
about 50 arcmin 1 2 in the GOODS-South field as a part of the WFC3 Early Release 
Science program. These LBGs at z ~ 1-3 are selected using dropout selection 
criteria similar to high redshift LBGs. The deep multi-band photometry in this 
field is used to identify best-fit SED models, from which we infer the following 
results: (1) the photometric redshift estimate of these dropout selected LBGs is 
accurate to within few percent; (2) the UV spectral slope /3 is redder than at high 
redshift ( z > 3), where LBGs are less dusty; (3) on average, LBGs at z ~ 1-3 
are massive, dustier and more highly star-forming, compared to LBGs at higher 
redshifts with similar luminosities, though their median values are similar within 
lcr uncertainties. This could imply that identical dropout selection technique, at 
all. redshifts, find physically similar galaxies; and (4) the stellar masses of these 
LBGs are directly proportional to their UV luminosities with a logarithmic slope 
of ~0.46, and star-formation rates are proportional to their stellar masses with 
a logarithmic slope of ~0.90. These relations hold true — within luminosities 
probed in this study — for LBGs from z ~ 1.5 to 5. The star-forming galaxies 
selected using other color-based techniques show similar correlations at z ~ 2, 
but to avoid any selection biases, and for direct comparison with LBGs at z > 3, 
a true Lyman break selection at z ~ 2 is essential. The future HST UV surveys, 
both wider and deeper, covering a large luminosity range are important to better 
understand LBG properties, and their evolution. 

Subject headings: galaxies: high redshift — galaxies: fundamental parameters — 
ultraviolet: galaxies — galaxies: evolution 


1. Introduction 

The high redshift frontier has moved to z > 7 as a result of the high resolution 
near-infrared (NIR) images from the Hubble Space Telescope (HST) Wide Field Camera 
3 (WFC3), and the Lyman break ‘dropout’ technique. The Lyman break technique was 
first applied to select Lyman break galaxies (LBGs) at z ~ 3 (Guhathakurta et al. 1990; 
Steidel etal. 1996, 1999), and since then it has been extensively used to select and study LBG 
candidates at redshifts z ~ 3-8 (e.g., Bouwens et al. 2007; Hathi et al. 2008b; Reddy & Steidel 
2009; Finkelstein etal. 2010; Yan etal. 2010). This dropout technique has generated large 
samples of faint star-forming galaxy candidates at z ~ 3-8. However, at highest redshifts 
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(z > 3) , it is very difficult to understand the details of their stellar populations using current 
space and ground-based telescopes. Their faint magnitudes make it extremely difficult to 
do spectroscopic studies, and limited high resolution rest- frame optical photometry make 
it challenging to investigate their spectral energy distributions (SEDs). These limitations 
make it imperative to identify and study LBGs at lower redshifts (z < 3). The primary 
reason for the lack of dropout selected LBGs at z ~ 1-3 is that we need highly sensitive 
space-based cameras to observe the mid- to near-ultraviolet (UV) wavelengths required to 
identify Lyman break at z ~ 1-3. The peak epoch of global star-formation rate at z ~ 1-3 
is now accessible using the dropout technique with the WFC3 UVIS channel. Hathi et al. 
(2010, hereafter H10) and Oesch etal. (2010) have used the HST WFC3 with its superior 
sensitivity to photometrically identify lower redshift (z ~ 1-3) LBGs. Understanding the 
LBGs at z < 3 is vital for two main reasons. First, we need to study the star- formation 
properties of these LBGs, because they are at redshifts corresponding to the peak epoch 
of the global star-formation rate (e.g., Ly etal. 2009; Bouwens etal. 2010). Second, they 
are likely lower redshift counterparts of the high redshift LBGs — because of their identical 
dropout selection and similar physical properties — whose understanding will help shed light 
on the process of reionization in the early universe (e.g., Labbe etal. 2010; Stark et al. 2010). 

There are primarily three techniques to select star-forming galaxies at z ~ 2: (1) sBzK 
(using the B , z, K bands, Daddi etal. 2004, 2007), (2) BX/BM (using the U, G, R bands, 
Steidel etal. 2004; Adelberger etal. 2004), and (3) LBG (using the bands which bracket the 
redshifted Lyman limit, H10; Oesch etal. 2010). All these approaches select star-forming 
galaxies, and yield insight into the star-forming properties of these galaxies, but they have 
differing selection biases, and so these samples don’t completely overlap (see Ly etal. 2011; 
Haberzettl etal. 2012, for details). Therefore, it is essential to apply identical selection crite- 
ria at all redshifts to properly compare galaxy samples and accurately trace their evolution. 
The LBG selection is widely used to select high redshift (z > 3) galaxies, and to do equal 
comparison with these galaxies, here we investigate physical properties of LBGs at z < 3. 

H10 used UV observations of the WFC3 Science Oversight Committee Early Release 
Science extragalactic program (PID: 11359, PI: O’Connell; hereafter “ERS”), which covers 
approximately 50 aremin 2 in the northern-most part of the Great Observatories Origins Deep 
Survey (GOODS; Giavalisco et al. 2004) South field, to identify LBGs at z ~ 1-3. The high 
sensitivity of the WFC3 UVIS channel data (Windhorst etal. 2011), along with existing 
deep optical data obtained with the Advanced Camera for Surveys (ACS) as part of the 
GOODS program are ideal to apply dropout technique in observed UV filters to select LBG 
candidates at z ~ 1-3. In this paper, we use this H10 sample of LBGs to investigate their 
physical properties by fitting stellar synthesis models to their observed SEDs. 
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This paper is organized as follows: In § 2, we summarize the WFC3 ERS observations, 
and discuss our LBG sample at 2 ~ 1-3 as well as the comparison sample of LBGs at 2 ~ 4-5. 
In 3 3, we fit observed SEDs of LBGs at 2 ~ 1-3 and 2 ~ 4-5 to stellar population synthesis 
models, and discuss the best-fit parameters (redshift, UV spectral slope, stellar mass, stellar 
age, and star- formation rates) obtained from these SED fits. In § 4, we discuss correlations 
between best-fit physical parameters and their implications on our understanding of LBGs. 
In § 5, we conclude with a summary of our results. 

In the remaining sections of this paper we refer to the HST /WFC3 F225W, F275W, 
F336W, F098M, F125W, F160W, filters as C/ 225 , U m , U 336 , Y 098 , J 125 , H m , to the HST /ACS 
F435W, F606W, F775W, F850LP filters as £ 435 , V 606 > ^ 775 , 2 850 , and to the Spitzer / IRAC 
3.6 pm, 4.5 pm filters as [3.6], [4.5], respectively, for convenience. We assume a Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) cosmology with G m =0.274, Ga = 0.726 and H 0 = 70.5 km 
s -1 Mpc -1 , in accord with the 5 year WMAP estimates of Komatsu etal. (2009). This cor- 
responds to a look-back time of 10.4 Gyr at 2 ~ 2. Magnitudes are given in the AB„ system 
(Oke & Gunn 1983). 


2. Observations and Sample Selection 

The WFC3 ERS observations (Windhorst etal. 2011) were done in both the UVIS (with 
a FOV of 7.30 arcmin 2 ) and the IR (with a FOV of 4.65 arcmin 2 ) channels. Here, we briefly 
summarize the UV imaging observations. The WFC3 ERS UV observations were carried 
out in three broad-band filters U 22 5 , U 275 and U 33 &. The C / 22 5 and C /275 filters were observed 
for 2 orbits (~5688s) per pointing, while the U 333 filter was observed for 1 orbit (~2778 
s) per pointing, for a total of 40 orbits over the full ERS field (8 pointings). We used the 
existing GOODS V 2 . 0 1 reduction of the ACS images in four optical bands (#435, Veso6, 2 775, 
2 8 5 o), which were re-binned to a pixel size of 0.09". To match the ERS IR (Y 098 , J125, //go) 
and re-pixellated ACS optical images, the UV mosaics have a pixel scale of 0.090" pix -1 
and cover ~50 arcmin 2 area of the GOODS-South field. Details of these observations and 
reduction process are described in Windhorst etal. (2011). 

The combination of the three WFC3 UV filters and the four ACS optical filters provide 
an excellent ability to select LBGs at 2 ~ 1-3 (H 10 ; Oesch etal. 2010), using the dropout 
technique to detect the Lyinan- break at rest- frame 912 A (Madau 1995). H 10 used dropout 
color selection technique in the ERS UV field to identify three sets of UV-dropouts — U225- 
dropouts, (/ 275 -dropouts and [/ 336 -dropouts — which are LBG candidates at 2 ~ 1.6, 2.2 and 


1 h.ttp:/ /archive.stsci.edu/pub/hlsp/goods/v2/ 



2 . 6 , respectively. They found 66 U 225 -, 151 U275- and 256 t/ 336 -dropouts to a magnitude limit 
of AB ~ 26.5 mag, respectively. 

In this paper, we start with the H10 LBG sample. Our goal is to investigate SEDs of 
reliable LBG candidates with at least 10-band HST coverage (augmented by additional data 
as described in § 3) from the H10 sample, so we apply the following filtering criteria. First 
criterion is the availability of the WFC3 IR (> 098 , J 125 , # 160 ) data. H10 used the WFC3 
UVIS and ACS data to select LBGs at z ~ 1-3. The WFC3 UVIS channel has a larger final 
ERS mosaic than the WFC3 IR channel, so we exclude LBG candidates that don’t have 
WFC3 IR data from our SED analysis. This criterion reduces the H10 sample size by about 
10%. Secondly, galaxies with poor SED fits (see § 3) as measured by their larger x 2 were 
excluded from the sample. This criterion removes additional ~10% of galaxies from the H10 
sample. This fraction of catastrophic x 2 outliers is consistent with the outlier fraction in 
the photometric redshift distribution of the dropouts in the H10 sample. The final sample 
of LBGs for the SED analysis — after applying above mentioned criteria — is 47 t/225-, 126 
t/ 2 75 - and 213 t/ 336 -dropouts. 

To compare SED properties of LBGs at z ~ 1-3, we select B 435 - and Veo6-dropouts in 
the WFC3 ERS field. These dropouts — LBG candidates at z ~ 3.7 and 4.7, respectively 
were selected following the Bouwens etal. (2007) selection criteria. The H10 sample is a 
comparatively brighter (~M* V ) sample, because of the limited depth and area of the ERS 
data, and hence, we also limit our comparison -B 435 - and V 606 -dropout sample to be around 
~M* V . Applying similar filtering criteria as LBGs at z ~ 1-3, we have a comparison sample 
of 155 R 4 35 -dropouts and 27 VWrdropouts in the ERS field. Based on Xue etal. (2011) 
X-ray catalog, there are four active galactic nuclei (AGN) in the H10 LBG sample (z < 3), 
three AGN in the B^-dropout sample, and none in the > 606 - dropout sample. These small 
number of X-ray AGN does not affect our results or conclusions. 

All subsequent analysis in this paper is done identically on these 5 samples ( U225- , t/ 275 -, 
t/ 336 -, B 435 -, and l^oe-dropouts) for proper comparison. To show general evolutionary trends, 
we combine three dropout samples from H 10 as a UV-dropout sample ( z — 1-3), and two 
high redshift samples as a B 435 -,Veo 6 -dropout sample (z ~ 4-5). 


3. Spectral Energy Distributions 

The Le PHARE software package (Arnouts et al. 1999; Ilbert et al. 2006) was used to mea- 
sure the photometric redshifts, and to fit the broadband SEDs of LBGs. The primary goal 
of SED fitting is to find the best-fitting synthetic stellar population model to the observed 
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photometry. Prom this best-fit model, we can estimate the redshift, stellar age, stellar mass, 
star-formation rate (SFR), dust extinction, and other physical properties of each galaxy. 
We use the 2007 version of the Bruzual & Chariot (2003, hereafter CB07) models, which 
has improved prescription of thermally pulsating AGB stars. We generated a set of stel- 
lar population models assuming a Salpeter initial mass function, and varying the redshift 
(z = 0. 1-6.0, <5z=0.1), metallicity (0.2, 0.4 and 1 Z 0 ), age (1 Myr < t < tu ), dust extinc- 
tion (0 < E(B-V) < 0.7 mag, using a modified Calzetti etal. 2000 attenuation law), and 
e-folding timescale (r=0.1-30 Gyr) for a star-formation history (SFH) ocexp(-t/r). The Le 
PHARE code assumes the Madau (1995) prescription to estimate inter-galactic medium (IGM) 
opacity. The model that gives the lowest x 2 is chosen as the best-fit SED. 

The contribution of major emission fines in different filters can be included in the models 
using the Le PHARE code. Neglecting emission fines during the SED fitting process can over- 
estimate the best-fit stellar ages and masses by about 0.3 dex (e.g., Schaerer Sz de Barros 
2009; Finkelstein etal. 2011; Atek etal. 2011). The Le PHARE code accounts for the con- 
tribution of emission lines with a simple recipe based on the Kennicutt (1998) relations 
between the SFR and UV luminosity, Ha and [Oil] fines. The code includes the Lya, Ha, 
H/9, [Oil], OIII[4959] and OIII[5007] fines with different line ratios with respect to [Oil] fine, 
as described in Ilbert etal. (2009). 

The observed photometry is available in up to 13 filters: three HST fWFCZ UVIS, four 
HST / ACS, three HST/ WFC3 IR, one VLT Ks, and two Spitzer/IRAC [3.6], [4.5] bands. 
We perform matched aperture photometry in 10 HST bands as discussed in H10, while 
we use VLT and Spitzer photometry from the publicly available GOODS-MUSIC catalog 
(Santini et al. 2009). The photometry in MUSIC catalog has accurate PSF-matching of space 
and ground-based images of different resolution and depth. Figure 1 shows example best-fit 
SEDs for LBGs at z ~ 1-3, and the comparison sample of LBGs at 2 ~ 4-5. 


3.1. Photometric Redshifts 

One of the free parameters during the SED fitting process is the redshift. To assess the 
accuracy of our SED-based photometric redshifts (z ph ) at z ~ 1-3, we compare them with 
the spectroscopic redshifts (z sp ) from various VLT/Magellan campaigns in the GOODS-S 
field (e.g., Grazian etal. 2006; Ravikumar etal. 2007; Vanzella etal. 2008; Wuyts etal. 2008; 
Balestra etal. 2010; Cooper etal. 2011). We find that only a small number (<30%) of our 
total dropout sample (H10) has spectroscopic redshifts, most likely due to the lack of strong 
features in 4500-9000 A range at 1 < z < 3, where most ground-based spectrographs on 
large telescopes are optimized. We matched 91 spectroscopic redshifts for the UV-dropout 
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sample (z ~ 1-3) selected with the criteria discussed in § 2. Figure 2 shows the comparison 
between the SED based photometric redshifts and the publicly available spectroscopic red- 
shifts. The catastrophic outhers — shown by concentric circles in Figure 2 — have quality 
flags that indicate the spectroscopic redshift is unreliable (in most catalogs C or worst). So 
it is likely that these spectroscopic redshifts are not correct and hence, redshift comparison 
for these objects is not credible. The histogram in Figure 2 shows the distribution of photo- 
metric redshift uncertainties 5z=(z ap — Zph/l+z sp ). Based on this distribution, we estimate 
<t(<Sz)~ 0.05, and < Sz >=-0.03. The fraction of catastrophic outliers (>3<r) is ~7%, exclud- 
ing objects with unreliable spectroscopic redshifts. Our photometric redshift uncertainties 
are consistent with Dahlen etal. (2010), who used the deepest and the most comprehensive 
photometric data in the GOODS-S field. Haberzettl etal. (2012) selected fairly bright LBGs 
at 2 ~ 2 using the GALEX data, and found similar photometric redshift uncertainties and 
ou;lier fraction for their dropout sample. The photometric redshift uncertainty in the im- 
plied redshift is also consistent with the dropout selection method applied to select these 
galaxies. The dropout selection technique uses the location of a spectral break within a 
photometric bandpass (filter), and therefore, the redshift uncertainty depends on the width 
of the bandpass, and could be as high as ~0.5 in z. 

The distribution of photometric redshift uncertainties (5z) is asymmetric, even after 
excluding objects with unreliable spectroscopic redshifts. There are more galaxies in the 
distribution with spectroscopic redshifts lower than their photometric redshifts i.e., (z ap — 
z p b < 0). Detail investigation of each ground-based spectra (if available) is needed to figure 
out what is causing spectroscopic redshift to be lower than photometric redshift. Such an 
investigation is beyond the scope of this paper, but we should point out that such asymmetric 
distribution is also observed for (TATEX-selected LBGs (e.g., Haberzettl etal. 2012), and is 
totally consistent within the estimated photometric redshift uncertainties. 


3.2. UV Spectral Slope 0 

The UV spectral slope 0 is determined from a power-law fit to the UV continuum 
spectrum (Calzetti etal. 1994), f\ oc A^, where f\ is the flux density per unit wavelength 
(ergs s -1 cm -2 A -1 ). We use the best-fit SEDs of dropout selected LBGs to estimate their 
UV spectral slope 0 by fitting a straight line between rest-frame 1300 and 1900 A in their 
model spectrum. This wavelength range covers 7 out of 10 spectral fitting windows identified 
by Calzetti etal. (1994) to estimate the UV spectral slope. This wavelength range is also 
ideal for comparing 0 values at higher redshifts, because those are usually measured between 
rest-frame 1600 and 2000 A. Figure 3 shows the slope-fitting method applied to the best-fit 
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SEDs to estimate 0, where the solid hue is the best-fit UV spectral slope, the dashed line 
is the best-fit SED, and the black filled circles are observed magnitudes. By selection (§ 2), 
we only consider galaxies with good SED fits so the choice of model should not affect the 
0 estimate. We estimate 0 for each galaxy, and then fit a Gaussian to the 0 distribution 
to find median (and sigma) value in each redshift bin. Table 1 shows median 0 values and 
their corresponding uncertainties for the UV-dropout, and the S 4 3 5 -,V 606 -dropout samples. 

The evolution in the UV spectral slope 0 as a function of redshift may indicate change 
in stellar populations of galaxies over cosmic time. We compare our 0 values with the higher 
redshift measurements from the literature (e.g., Bouwens etal. 2011; Finkelstein etal. 2012). 
Figure 4 shows the UV spectral slope 0 as a function of redshift. Blue filled squares are 
median 0 values measured between rest-frame 1300 and 1900 A for our dropout samples. 
To test how 0 measurements are affected by the selection of rest-frame UV wavelength 
range, we also measured 0 between rest-frame 1300 and 3400 A, which are shown by blue 
open squares in Figure 4. Both 0 values agree within la uncertainties. The red filled 
diamonds are measurements from Bouwens etal. (2009, 2011) and purple filled circles are 
from Finkelstein etal. (2012). The 0 values from Bouwens etal. (2009, 2011) are for the 
galaxies around M* v , which is consistent with our sample, while Finkelstein etal. (2012) 
measurements are based on all galaxies extending to those fainter than M* v in their respective 
redshift bins. The uncertainties on the median values of 0 are the standard error of the mean 
in the case of Finkelstein etal. (2012) and our measurements, while Bouwens etal. (2009, 
2011) uncertainties represent ler scatter. For comparison, our estimated ler scatter in median 
0 values are listed in Table 1. 

Figure 4 shows that the median values of 0 decreases as redshift increases (0 ~-1.6 
at z ~ 1.6 to 0 ~-2.4 at z ~ 8), which could imply variations in one or more physical 
properties of LBGs as a function of redshift. Figure 5 shows evolution in 0 as a function 
of best-fit SED parameters (stellar mass, stellar age, dust content, SFR) and redshift. The 
lowest redshift bin {z ~ 1.6) is shown by the smallest circles, and the highest redshift bin 
(z ~ 4-5) is shown by the largest circles. The largest change (factor of ~2 or 0.3 dex) is 
seen in the dust content E(B- V) of galaxies as 0 changes from -1.6 (at z ~ 1.6) to -1.9 (at 
2 ~ 4-5), while other three parameters vary much less than a factor of 2. This could imply 
that change in the dust content of LBGs has largest effect on the UV spectral slope 0 , and 
any variation in 0 as a function of redshift could most likely be due to changing dust content 
of galaxies. Therefore, based on our 0 estimates, as shown in Figure 4, we could say that 
LBGs at lower redshift ( z ~ 1.6) have more dust than LBGs at higher redshift (z ~ 5). This 
trend of 0 is consistent with those from previous studies, which have found that galaxies at 
z ci : 6 tend to be bluer than those at z ~ 3 (e.g., Stanway etal. 2005; Bouwens etal. 2006; 
Hathi etal. 2008a; Wilkins etal. 2011). Those 0 measurements, on uniformly selected LBGs, 
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were limited to LBGs at 2 > 3, and our results below z~3 extends this observed trend to 
z ~ 1.5. 

The evolution in the UV spectral slope j3 could also be due to changing star-formation 
history, initial mass function (IMF), and/or metallicity. These effects are believed to be 
much smaller than the effects from changing dust content of the galaxy. Many authors 
have investigated various stellar population models to estimate these effects. Bouwens et al. 
(2011) explored sensitivity of the UV-continuum slope j3 to changes in the mean metallicity, 
age, or dust extinction by choosing one fiducial model as a benchmark, and then changing 
various model parameters to assess changes in /?. They conclude that a factor of 2 (or 0.3 dex) 
changes in metallicity, age or dust content result in 0.07, 0.15, 0.35 changes in the UV spectral 
slope /3, respectively. This implies that changes in the dust content have much larger effect 
on the UV-continuum slope than similarly-sized changes in the age, metallicity, or the stellar 
IMF. Similar studies (e.g., Leitherer etal. 1999; Hathi etal. 2008a; Wilkins etal. 2011) have 
come to the same conclusion that though (3 could be affected by various stellar population 
properties, the change in dust content of galaxies is the predominant effect which causes /3 
to change. We should also emphasize that it is very challenging to completely understand 
these various effects based on observations only, rigorous modeling and/or simulations are 
required to fully assess the contributions of these various effects on the UV spectral slope ft. 


3.3. Stellar Population Properties 

We use the 10-band HST photometry along with the VLT Ks and the Spitzer /1RAC 
[3.6], [4.5] bands to construct an observed SED which we then compare to a suite of model 
templates from CB07 and find the best-fit model through y 2 minimization. The model 
parameters are outlined in § 3. The best-fit model allows us to estimate stellar age, stellar 
mass, dust extinction E(B-V) and SFRs for each galaxy. Uncertainties (~0.3-0.4 dex) 
in stellar ages, masses and SFRs are estimated based on the uncertainties in redshift and 
observed photometry. 

One of the main limitations of SED fitting is the need to assume a SFH, which cannot 
be reliably constrained for a galaxy from limited photometric data points. We have assumed 
an exponentially declining SFH. Different SFHs (e.g., rising, constant, declining) introduce 
systematic uncertainties in the stellar mass determinations, mostly at redshift greater than 
z ~ 3 (e.g., Lee etal. 2010; Papovich etal. 2011). These uncertainties are typically <0.3 
dex (e.g., Finlator etal. 2007), and are within our estimated uncertainties. Stellar ages are 
highly sensitive to the assumed SFH. Any prior star- forming episode can be overshadowed by 
newly born stars from the most recent star-formation, totally neglecting possible existence 



of older population in a given galaxy. Therefore, based on assumed histories, it is possible 
to get an older or a younger age for the same galaxy. Hence, interpreting the stellar ages 
derived from SED fitting can be tricky, and at the very least, the uncertainties on the stellar 
ages could be much larger than estimated uncertainties (~0.3-0.4 dex). Because of these 
issues, in subsequent analysis, we will not elaborate on stellar population ages and focus on 
other physical properties. 

Figure 6 shows the distributions of stellar age, stellar mass, SFRs, and E(B-V) for 
LBGs at z ~ 1-3 (black), and the comparison sample of LBGs at z ~ 4-5 (red). The 
median values are shown by dashed vertical lines and la uncertainties in these distributions 
for LBGs at z ~ 1-3 are shown by an error bar at the top of the black histogram. A two- 
sided K-S test indicates a probability less than 0.006 that the distributions (red and black 
histograms) are drawn from the same parent distribution. Figure 6 shows a general trend 
that — on average — higher redshift LBGs have low SFRs, less dust, and are less massive 
than their lower redshift counterparts, though median values of two distributions (red and 
black) are similar within la uncertainties. This result is in good agreement with previous 
studies comparing LBGs at .z ~ 3 and z ~ 5 (e.g., Verma etal. 2007). The average E(B-V) 
at z ~ 2 is consistent with studies based on star-forming galaxies selected using BX/BM 
color technique (e.g., Erb etal. 2006; Sawicki 2012). The distribution of E(B-V) completely 
agrees with the UV spectral slope evolution as discussed in § 3.2, implying that the LBGs 
at z ~ 1-3 are more dusty (redder) compared to LBGs at z~ 4-5. 


4. Results and Discussion 
4.1. Stellar Mass vs UV Luminosity Relation 

The rest-frame UV light traces current or instantaneous SFR, while rest-frame optical 
and NIR data help us to estimate stellar masses of galaxies. If the galaxy stellar mass and 
UV luminosity are related then we can directly use rest-frame UV light to estimate stellar 
mass without needing rest-frame optical/NIR data. Figure 7 shows stellar mass of LBGs at 
2 ~ 1.5-5 as a function of their UV absolute magnitude. These quantities are based on best- 
fit 3EDs, and their typical uncertainties are shown in the lower-left corner. The dotted lines 
are best-fit line obtained by keeping the logarithmic slope fixed at 0.46, which was estimated 
by Sawicki (2012) for star-forming galaxies at z ~ 2. The dot-dash lines show the scatter 
from the best-fit line, which is ~0.3 dex for LBGs at z ~ 1-3 and about 0.2 dex for LBGs 
at z ~ 4-5. We also tested the validity of this relation by fitting the slope of the line rather 
than fixing it. We find that the fitted slope is in the range of 0.42±0.06 for our LBG samples, 
which is consistent with 0.46 within the estimated la scatter in this relation. Therefore, we 
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find that a proportionality relation between these two parameters with a logarithmic slope 
of 0.46 provides a good fit to the data. The stellar masses of the brighter LBGs — with 
UV luminosities near the L* v value of LBGs at z ~ 3 from Steidel etal. (1999) — are about 
a factor of 2 lower than 10 10 M e estimated by Papovich etal. (2001). This discrepancy, 
though within our estimated uncertainties, could be due to the fact that we include emission 
lines in our SED fitting which can affect stellar masses by a factor of ~2. The stellar mass- 
UV luminosity relation is fairly tight with a small scatter (<0.3 dex), which is consistent 
with other studies at similar redshifts (e.g., Papovich etal. 2001), and it points to a nearly 
constant mass-to-light ratio (log(M/L) ~-0.5, see § 4.3) for LBGs between z ~ 1.5 and 5. A 
similar correlation between stellar mass and absolute magnitude has been reported for LBGs 
at z ~ 5-6 by Stark etal. (2009). 

Figure 7 shows that LBGs at z ~ 1.5-5 follow similar linear correlation between stellar 
mass and UV absolute magnitude (within uncertainties) for M uv between -19 and -22.5 mag. 
It is important to note that Shapley et al. (2005) does not find any correlation between the 
stellar mass and UV absolute magnitude for star-forming galaxies at z ~ 2 with stellar 
masses > 10 10 M®. This could be due to different color-selection technique (BX/BM) used 
by the Shapley et al. (2005) to select star-forming galaxies, whose physical properties could 
differ from the dropout selected LBGs at these masses (e.g., Ly etal. 2011; Haberzettl etal. 
2012). It is also possible that their sample, which consists of spectroscopically confirmed 
bright galaxies with stellar masses greater than or equal to 10 10 M®, has more massive 
galaxies than our sample because of the limited survey area in the WFC3 ERS field. The 
ERS observations are too limited in area and depth to cover a larger luminosity range, so 
we cannot predict how this relation will evolve for luminous (M uv < -22.5 mag) or dwarf 
(A/ uv > -19 mag) galaxies at these redshifts. 


4.2. SFR vs Stellar Mass 

The correlation between the current SFR and stellar mass in star-forming galaxies, 
also known as ‘main sequence of star-formation’ (MS), has been observed at z < 2 (e.g., 
Noeske etal. 2007; Elbaz etal. 2007; Daddi etal. 2007). These studies have shown that the 
MS relation seems to be not evolving strongly with redshift, but the zeropoint does: that is 
high redshift (z — 2) star-forming galaxies arc forming stars at a higher rate than similar 
mass local galaxies. In Figure 8, we investigate this relation and star- formation histories 
for LBGs at z ~ 1.5-5. These quantities are based on best-fit SEDs, and their typical 
uncertainties are shown in the lower-right corner. The dotted lines are best-fit line obtained 
by keeping the logarithmic slope fixed at 0.90, estimated for star-forming galaxies at z < 2 



(e.g., Elbaz etal. 2007; Daddi etal. 2007; Sawicki 2012). The dot-dash lines show the scatter 
from the best- fit line, which is ~0.6 dex for LBGs at z ~ 1-3 and about 0.4 dex for LBGs 
at z ~ 4-5. We also obtained the best-fit logarithmic slope for this relation, and found 
the fitted slope in the range of 0.81±0.30 for our LBG samples, which is consistent with 
0.90 within the estimated la scatter in this relation. We find that a proportionality with 
a logarithmic slope of 0.90 provides a good fit to the data with few outliers at stellar mass 
greater than 10 10 M©. 

Finlator et al. (2006) have shown that tight relation exists between SFR and stellar 
mass for galaxies at z ~ 4 using the cosmological hydrodynamic simulations, which is also 
consistent with the observations (Bouwens etal. 2011). Finlator etal. (2006) also point out 
that the scatter in the Figure 8 could be a measure of SFR ‘burstiness’ as a function of stellar 
mass. This means that the linear relation (with a logarithmic slope of ~0.90) indicate an 
average SFR for a given stellar mass, but galaxies can also experience bursts of up to two 
times the average SFR value at the same stellar mass as shown by the scatter. The scatter 
in the SFR versus stellar mass relation for LBGs at z ~ 1-3 is slightly larger than ~0.3 dex 
— observed at z ~ 2 by Daddi et al. (2007) — possibly because of few galaxies forming a 
sharp edge towards high SFR values, as shown in the relation for U 225 - and E/ 275 -dropouts 
(upper panel in Figure 8). These galaxies have low stellar ages (less than 10 Myr), which 
could be highly uncertain as discussed in § 3.3. It is also possible that this edge could be an 
artifact due to lower limits on the model parameters r and t (e.g., Hainline et al. 2012). We 
also note that McLure etal. (2011) argue that the tightness in the SFR-stellar mass relation 
depends on the assumed SFH. The scatter in this relation is much less for a constant SFH, 
while it is much larger for other SFHs. Therefore, it is also likely that the larger scatter we 
see in Figure 8 could be due to different SFHs. 

Figure 8 shows that, though our data has little more scatter compared to the MS 
relation at z < 2, the majority of our galaxies fall on to this relation characterized by a 
logarithmic slope of 0.90. A similar correlation is observed at z cz 6-8 by McLure et al. 
(2011), and supported by cosmological hydrodynamic simulations of Finlator etal. (2011). 
Our observations confirm this MS relation for star-forming galaxies from z ~ 1.5 to 5, 
implying that — on average their star-formation histories are similar. 


4.3. Implications 

In previous sections, we have shown that LBGs at z ~ 1-3 — on average — are massive, 
dustier, and have higher star-formation rates than LBGs at z~ 4-5 with similar luminosi- 
ties. though it should also be noted that they are not very different within estimated la 
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uncertainties. As pointed out by Papovich etal. (2011), the number densities of galaxies 
at fixed luminosity could change substantially over this redshift range, which could lead to 
potential biases when comparing galaxies at different redshifts. However, the general trends 
we observe in stellar masses, SFRs, and dust extinction are supported by other independent 
means. The characteristics UV luminosity (L* v ) is increasing as a function of redshift from 
z ~ 8 to 2 (e.g., H10), which implies increase in SFRs with time, while Finkelstein etal. 
(2010) have shown that stellar masses for M* v LBGs grow from 2 ~ 8 to 2. The UV spectral 
slope (3 shows evolution as a function of redshift (Figure 4) , which could indicate lower dust 
content at higher redshifts. The higher dust content in LBGs at lower redshift is also in 
accordance with the studies at z ~ 1 (e.g., Burgarella et al. 2007; Basu-Zych etal. 2011), 
while the Verma et al. (2007) supports the lower dust content in LBGs at 2 ~ 5. Therefore, 
the ensemble properties of LBGs in our sample are in general agreement with the expected 
results. 

The stellar mass-UV luminosity relation (Figure 7) and the SFR-stellar mass relation 
(Figure 8) are based on measurements from best-fit SEDs, therefore, it is possible that these 
quantities are not totally independent which might affect their observed correlations. To 
investigate this, we show distributions of mass-to- light (M/L; Mass/L uv ) ratios and specific 
SFRs (SSFR; SFR/Mass) in Figure 9. The black (red) histograms show distribution for 
LBGs at 2 ~ 1-3 (2 ~ 4-5), and the median values are shown by dashed vertical fines. The 
median values of M/L ratio and SSFR for LBGs at 2 ~ 4-5 are slightly lower than that at 
2 ~ 1-3, but are still consistent within the la uncertainties as shown by the error bar on the 
top of the black histogram. A two-sided K-S test indicates a probability less than 0.05 that 
the distributions (red and black histograms) are drawn from the same parent distribution. 
The constancy of the M/L ratio and SSFR between 2 ~ 1.5 and 5 agrees very well with the 
constant slope we find in Figure 7 and Figure 8 for our sample of LBGs. 

Stellar masses of LBGs at 2 ~ 1-3 are generally well correlated with UV absolute 
magnitude and current SFR, as expected for star-forming galaxies at similar redshifts (e.g., 
Elbaz etal. 2007; Daddi etal. 2007; Sawicki 2012). These correlations implies very similar 
mass assembly and SFH for these galaxies, but the exact nature of SFHs is still not clearly 
understood. Papovich etal. (2011) showed that the cosmologically averaged SFRs of star- 
forming galaxies at 3 < 2 < 8 — at constant co-moving number density — increase smoothly 
from 2 = 8 to 3, and the stellar mass growth in these galaxies is consistent with this derived 
SFH. The scenario of rising SFH (see also Lee et al. 2010) is also supported by recent results 
from the cosmological hydrodynamic simulations (e.g., Finlator et al. 2011). The models with 
rising SFHs conflicts with the assumptions that the SFR in distant galaxies is either con- 
stant or decreasing exponentially with time (e.g., Papovich etal. 2005; Shapley etal. 2005; 
Labbe etal. 2010). Though, we remind the reader that the models with rising SFHs advo- 
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cated by Papovich etal. (2011) and others correspond to a cosmologically averaged SFHs 
for typical galaxies, and not individual galaxies, because they could involve random events 
that changes their instantaneous SFR. Papovich etal. (2011) and Lee etal. (2010) also argue 
that rising SFHs are most beneficial to higher redshift ( z > 3) galaxies. We find that for our 
assumed SED model parameters, the LBGs between redshift z ~ 1.5 and 5 — on average — 
have similar SFHs, though the precise nature of SFHs at all redshift is still under debate, 
and could also affect the SFR- stellar mass correlation. 

Our analysis demonstrates that the dropout selected galaxies at z ~ 1-3 — within lumi- 
nosities probed here — show similar correlations between physical parameters (SFR, stellar 
mass, UV luminosity) as other star-forming galaxies selected using different color criteria 
(e:g., sBzK , BX/BM) at z ~ 2. This is consistent with the Lyetal. (2011) conclusion 
that majority (~80 -90%) of the dropout selected galaxies overlap with other color selected 
star-forming galaxies with stellar masses less than 10 5 * * * * 10 M©. The stellar mass range for our 
current sample is between ~10 8 and ~10 10 M 0 , with sample completeness around 10 9 9 5 

M(t. Significant differences between the dropout selected sample and other color selected 
samples of star-forming galaxies at z ~ 2 exists for massive galaxies (> 10 10 M 0 ; Ly etal. 
2011). Therefore, it is vital to use uniform selection technique at all redshifts to avoid any 
selection biases. The Lyman break dropout technique is the most convenient and widely 
used method to select galaxies at z > 3, and we have shown that LBGs at z ~ 1-3 selected 
using this dropout technique have similar physical properties (within uncertainties) as LBGs 
at z ~ 4-5 with similar luminosities. Hence, LBG selection at z < 3 is important to under- 
stand properties of LBGs and properly investigate their evolution as a function of redshift. 
The validity of LBG properties over wide luminosity and mass range can be investigated in 
detail with the upcoming and future WFC3 UV surveys such as CANDELS (Grogin et al. 
2011; Koekemoer etal. 2011) and the WFC3 UV UDF (Rafelski etal. 2012, PI: H. Teplitz). 

5. Summary 

In this paper, we investigated stellar populations of LBGs at z ~ 1-3 selected using 

HST/ WFC3 UVIS filters in the GOODS-S field. We used deep multi- wavelength observa- 
tions from the HST, VLT, and Spitzer to compare observed SEDs with the spectral synthesis 

models to infer physical properties (stellar masses, stellar ages, SFRs, and dust extinction) of 
these LBGs. We also compared these LBGs with their higher redshift (z ~ 4-5) counterparts. 

Our results can be summarized as follows: 

• We obtain reliable photometric redshifts for dropout selected LBGs at z ~ 1-3 based 
on 10-13 band SEDs. 



• The UV continuum slope j3 for LBGs at z ~ 1-3 is redder compared to their higher 
redshift counterparts, implying higher dust content in these LBGs. 

• On average, LBGs at z ~ 1-3 are massive, dustier and more highly star-forming com- 
pared to LBGs at z ~ 4-5, though their median values are very similar within estimated lcr 
uncertainties. This similarity emphasizes the importance of identical Lyman break selection 
technique at all redshifts, which selects physically similar galaxies. 

• The stellar mass-absolute UV magnitude relation for LBGs between z ~ 1.5 and 
5 show linear correlation with a logarithmic slope of ~0.46, while the SFR-stellar mass 
relation show similar correlation with a logarithmic slope of ~0.90. To properly compare 
and interpret such relations at higher {z > 3) redshift, and to avoid any selection biases due 
to different selection techniques, a true Lyman break selection is required at 2 ~ 2. 

• We need larger HST UV surveys to cover full range in luminosity/mass and better 
understand LBG properties, and their evolution. Both deeper and wider UV surveys are 
needed. The wider one to probe the high mass end, while the deeper one will probe the sub- 
L* population. A large number of HST orbits have been used for dropout selected galaxies 
at 2 > 3, and the lower redshift regime needs to be explored in a comparable manner. 


This paper is based on Early Release Science observations made by the WFC3 Scien- 
tific Oversight Committee. We are grateful to the Director of the Space Telescope Science 
Institute for awarding Director’s Discretionary time for this program. Support for program 
#11359 was provided by NASA through a grant HST- GO- 11 359. 08- A from the Space Tele- 
scope Science Institute, which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. This research was (partially) sup- 
ported by a grant from the American Astronomical Society. 
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Fig. 1. — [Top panel] Example best-fit SEDs of LBGs at z ~ 1-3. We have used the Le 
PHARE code (Arnouts etal. 1999; Ilbert etal. 2006) to compute photometric redshifts and to 
perform SED fitting. Grey concentric circles are observed magnitudes in three HST/ WFC3 
UVIS, four HST /AGS, three HST/ WFC3 IR, one VLT Ks, and two Spitzer/TRAC [3.6], 
[4.5] bands. [Bottom panel] Same as the top panel but for the comparison sample of LBGs 
at z ~ 4-5. Stellar masses are in solar units and stellar ages are in years. 
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Fig. 2. — Comparison of the SED-based photometric redshifts with the publicly available 
spectroscopic redshifts for LBGs at 2 ~ 1-3. The concentric circles indicate unreliable 
spectroscopic redshifts based on their quality flags. The dotted line shows the 1-to-l relation 
between two redshifts. The histogram shows distribution of ( z sp — z p h/l+z sp ), and the la 
uncertainties in this distribution is ~5%. The dot-dash lines show 3 <7 limits. 
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Fig. 3. — Example fitting of the UV spectral slope to best-fit SEDs for two LBGs from the 
2 ~ 1-3 sample. The dotted black line is the best-fit SED, while solid black line shows the 
estimated UV continuum slope fi. The black filled circles are observed magnitudes. We test 
two different rest-frame UV wavelength range (1300-1900 A and 1300-3400 A) to address 
the robustness of our estimated UV continuum slope /3, and find similar median slopes within 
la scatter. 
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Redshift 

Fig. 4. — UV spectral slope /3 as a function of redshift. Blue filled squares indicate the 
median f} values when we fit rest-frame 1300 to 1900 A wavelengths, arid blue open squares 
indicate /?s when we fit 1300 to 3400 A. The higher redshift (z > 2.5) measurements from 
Finkelstein etal. (2012) and Bouwens etal. (2009, 2011) are also shown for comparison. Ours 
and the Finkelstein etal. (2012) uncertainties are the standard error of the mean, while 
uncertainties from Bouwens et al. are la scatter. For comparison, our la uncertainties are 
listed in Table 1. 
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Fig. 5. — UV spectral slope fl versus best-fit SED parameters (E(B-V), stellar mass, stellar 
age, star-formation rate). The size of the circle increases with redshift i.e., the smallest circle 
correspond to the lowest redshift bin ( z ~ 1.6), and the largest circle correspond to z ~ 4-5 
bin. The changes in the UV-continuum slope (3 from z ~ 1.6 to 5 are most likely due to 
the change (factor of ~2 or 0.3 dex) in the dust content E(B-V) of the galaxies, as other 
parameters show smaller variation as a function of redshift. 





Fig. 6. — Distribution of stellar masses, stellar ages, star-formation rates, and dust extinction 
for LBGs at z ~ 1-3 (black) and the comparison sample at z ~ 4-5 (red). All four parameters 
are estimated from the best-fit SEDs. Median values axe shown by vertical dashed lines and 
1 <t uncertainties are shown by an error bar at the top of the black histogram. Though the 
higher redshift LBGs have, on average, lower values of most of these physical parameters 
than LBGs at z ~ 1-3, the median values are similar within estimated uncertainties. Stellar 
ages do not show any clear evolution with redshift mainly because of higher uncertainties in 
their measurements (see § 3.3). The histograms show normalized numbers for both samples 
(total number of galaxies in each sample is shown in Table 1), and a two-sided K-S test 
indicates a probability P < 0.006 that the two distributions are drawn from the same parent 
distribution. 
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Fig. 7. — Stellar mass versus UV absolute magnitude relation for 7 / 225 -, 7/275-, 7 / 336 -dropout 
samples. Bottom right panel shows the same relation for the comparison sample of LBGs at 
z ~ 4-5 (black circles for z ~ 4, and black squares for z ~ 5). The average uncertainties are 
shown in the bottom-left corner. The dotted black line is the best-fit line with a logarithmic 
slope of 0.46. The dot-dash, line shows the lcr scatter (~0.3 dex) from the best-fit linear 
relation. 
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Fig. 8. — Stellar mass versus SFR relation for U 225 -, U 275 -, [/336-dropout samples. Bottom 
right panel shows the same relation for the comparison sample of LBGs at z ~ 4-5 (black 
circles for z ~ 4, and black squares for z ~ 5). The average uncertainties are shown in the 
bottom-right corner. The dotted black line is the best-fit line with a logarithmic slope of 
0.90. The dot-dash line shows the la scatter (~0.6 dex, ~0.4 dex for z ~ 4-5) from the 
best-fit linear relation. In the upper panel, few galaxies form a sharp edge towards high SFR 
values, which could be an artifact due to lower limits on the model parameters r and t (e.g., 
Hamline etal 2012). 
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Fig. 9. — Distribution of mass-to-light ratios and specific SFRs for LBGs at 2 ~ 1-3 (black) 
and for the comparison sample at z ~ 4-5 (red). Median values are shown by vertical dashed 
lines and lcr uncertainties are shown by an error bar at the top of the black histogram. 
Though the higher redshift LBGs have, on average, lower values of these parameters than 
LBGs at ^ ~ 1-3, the median values are similar within estimated uncertainties. The his- 
tograms show normalized numbers for both samples (total number of galaxies in each sample 
is shown in Table 1), and a two-sided K-S test indicates a probability P < 0.05 that the two 
distributions are drawn from the same parent distribution. 
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Table 1. UV Spectral Slopes (/?) 


Redshift 
< z > 

Number of 
Galaxies (N) 

UV Slope 0 

P 

1 a 

scatter 

SEM 6 

1.6 (±0.2) 

47 

-1.59 (-1.67) 

0.29 (0.27) 

0.04 (0.04) 

2.2 (±0.3) 

126 

-1.71 (-1.81) 

0.34 (0.37) 

0.03 (0.03) 

2.6 (±0.4) 

213 

-1.71 (—1.78) 

0;47 (0.50) 

0.03 (0.03) 

3.8 C (±0.5) 

182 

-1.88 (-1.93) 

0.45 (0.41) 

0.03 (0.03) 


a The UV continuum slope /3 is estimated from best-fit SEDs by fitting 
a line between the rest-frame wavelengths 1300-1900 A (1300-3400 A) 

b The standard error of the mean = (lcr scatter) /y/N 

c The comparison sample of B 435 - and V 606 -dropouts. Because of small 
number of Vfioe-dropouts, the average/median redshift is similar to a 
^435-dropout. 





